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PREFACE 

The  main  purpose  of  Tn  No.  64T-1  is  to  serve  as  a  guide  to  students  ■ 
using  the  Fortran  cards  for  computer  program  "CMPRSOR"  to  reduce  the 
data  taken  during  a  classroom  test  of  the  compressor  test  rig.   This 
compressor  test  rig  is  described  in  USNPGS  Technical  Report  No.  12. 

To  reduce  data  taken  during  a  compressor  test  rig  run,  it  is 
necessary  only  to  follow  directions  on  pages  7S  8  and  9  of  this  technical 
note,   Appendix  A  will  also  help  toward  an  understanding  of  the  program 
and  its  use. 

Appendix  B  gives  the  information  necessary  for  slide  rule  calculations 
of  efficiency  data. 

In  content,  TN  No.  64T-1  is  a  compilation  of  information  revelant 
to  the  compressor  test  rig  not  contained  in  USNPGS  Technical  Report  No.  12 
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SYMBOLS : 

A 

Amps. 

S 

Slip  -  RPM 

TF 

Temperature  -  ©F 

Tc 

Temperature  =  C 

V 

Volts 

W 

Watts 

Wl 

Watts  (Left  Hand  Watt  meter) 

W2 

Watts  (Right  Hand  Watt  meter 

wT 

Watts  Loss 
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2.  INTRODUCTION: 

Technical  Report  No.  12  published  in  February  1955  by  Professors 
Vavra  and  Gawain  covers  the  construction  and  use  of  the  compressor 
test  rig  and  gives  all  information  required  for  data  reduction  leading 
to  values  for  flow  and  pressure  coefficients  cp  and  f.   The  derivation 
fo  formulas  for  flow  and  pressure  coefficients  is  also  included  in  TR  No.  12. 
The  information  given  for  slide-rule  computation  of  flow  and  pressure  co- 
efficients cp  and  \|f  is  still  pertinent;  however,  the  computer  program  des- 
cribed and  discussed  in  the  above  report  applied  to  a  computer  no  longer 
at  the  USNPGS. 

The  following  supplements  Report  No   12,  in  two  respects :   First s 
it  shows  the  development  of  formulas  for  calculation  of  compressor 
efficiency  by  the  use  of  a  computer  program  andj second,  it  gives  detailed 
information  on  the  use  of  a  Fortran  program.   This  program  is  for  the 
reduction  of  all  data  taken  during  a  compressor  test  rig  run. 

3.  COMPUTER  PROGRAMMING  OF  MOTOR  OUTPUT  DATA: 

Before  1956  the  compressor  test  rig  was  operated  without  power  input 
instrumentation.   In  1956  the  Westinghouse  Company  calibrated  the  compressor 
drive  motor  and  supplied  loss  curves. 

Net  input  to  the  compressor  is  equal  to  the  input  to  the  motor  minus 
losses,  since  the  compressor  rotor  is  connected  directly  to  the  motor  and 
the  motor  bearings  carry  all  bearing  loads.   Input  to  the  motor  is  the 
algebraic  (since  one  wattmeter  may  have  a  negative  reading)  sum  of  readings 
from  two  wattmeters,  the "two -wattmeter  system"  for  3-phase  current.   Total 
losses  may  be  compiled  by  using  the  loss  curves.   Compressor  efficiency, 
then,  is  the  theoretical  power  required  by  the  compressor  divided  by  the 
net  input  to  the  compressor.   A  sample  calculation  (for  slide-rule  computa- 
tion) of  the  motor  output  and  curves  from  data  supplied  by  the  Westinghouse 
Company  are  given  in  Appendix  B 

In  order  to  use  the  curves  as  part  of  a  computer  program  it  was  necessary 

to  reduce  them  to  formulas   Two  of  the  curves  are  straight  lines.   Another 

curve  has  only  three  points.   It  was  assumed  that  this  curve  is  parabolic 

so  the  formula  for  the  curve  was  set  up  as: 

2 
y  -  A  +  BX  +  CX 
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or,  specifically: 

2 
Watts  Loss  =  A  +  BX  Volts  +  CX  (Volts) 

Substituting  values  at  the  three  points  given,  three  equations  can  be 
written  from  which  the  unknowns  A,  B,  and  C  can  be  computed. 

The  remaining  curve  has  seven  data  points.   With  the  exception  of 
one  point,  these  points  are  on  a  parabolic  curve  determined  by  any 
three  of  the  six  points.   The  points  were  checked  and  found  to  lie 
on  this  theoretical  curve  to  an  accuracy  greater  than  the  accuracy 
to  which  the  curve  drawn  through  the  points  can  be  read. 

For  values  above  12  amperes  the  constants  of  another  parabolic 
curve  were  computed,  using  the  12  ampere  point,  the  14  ampere  point  and 
a  point  obtained  by  graphical  interpolation  between  the  12  ampere  point 
and  the  14  ampere  point. 

The  computer  program  calls  for  the  computer  to  choose  which  formula 
to  use,  depending  on  whether  the  input  is  above  or  below  12  amperes-   With 
Fortran  a  choice  is  called  for  by  the  "if"  statement. 

A  program  for  motor  input  and  loss,  and  thus  for  compressor  efficiency, 
was  first  written  to  supplement  the  original  compressor  test  rig  program 
written  for  the  CDC  102.   Reduction  of  efficiency  data  is  now  part  of 
a  complete  data  reduction  program  for  the  compressor  written  in  Fortran. 

For  slide  rule  computations: 

Efficiency  ■   Theoretical  Power  of  Compressor/Input  to  Compressor 

Flow  x  Pressure  x  A  Constant/Electrical  Input  -  (Motor  Losses  + 

TR) 

"T  "  is  the  power  loss  due  to  the  resisting  torque  of  the  water  seal 
R 

mechanism. 

Input    to    the  motor    is  measured   by   two   wattmeters 

W   =  W]_   +  W^ 

At  low  power  input  the  right  hand  wattmeter  (W~)  may  be  negative  ; 
if  this  is  indicated,  two  leads  are   reversed  and  the  reading  is  recorded 
as  negative. 

For  slide  rule  computations  the  meter  readings  must  be  multiplied 
by  a  factor  to  give  true  volt,  ampere,  and  watt  values.   For  use  in 
compiling  data  fotr   the  Fortran  program  "CMPRSOR"   meter  readings  are 
recorded  exactly  as  read. 
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The  ampere  meter  has  multiple  scales.   The  correct  amperage  is 
obtained  by  reading  the  second  scale  from  the  bottom  and  multiplying  by 
10.   For  slide  rule  computations  an  average  istaken  of  the  readings 
from  the  three  meters. 

Multiplying  Factors 

Amperes  (Second  Scale  from  Bottom)   Multiply  by  10 

Volts  n    "   20 

Watts  "    "   40 

For  compiling  motor  losses  reference  may  be  made  to  the  loss  curves 
and  sample  calculations  in  Appendix  A. 

As  an  alternate  method  of  computing  motor  losses  for  slide  rule 
calculation  of  compressor  efficiency  the  loss  curve  formulas  may  be 
used. 

The  combined  friction,  windage  stray  and  core  loss  is  a  function 
of  the  voltage. 

WI    =     -210*.  0  +  0.8929   V  -   0.000104/  V2 

Stator  copper  losses  are  a  function  of  the  amperage. 

WT  =  2.73  A2  when  A  <  12 

W  =  4484  -  783A  +  33. 0A2 

La 

when  A  >  12 

This  loss  is  modified  by  a  factor  which  is  a  function  of  motor  winding 
temperature. 

Factor  =  0.904  +  0.00385  Tc 
Rotor  copper  losses  are  an  indirect  function  of  the  slip.   To  get 
these  losses  the  gross  watts  input  minus  the  stator  losses  is  multiplied 
by  a  factor  which  is  a  function  of  slip. 

Factor  =  0.000556  x  Slip 

When  the  compressor  is  set  up  for  taking  rotor  blade  pressures  there 
is  an  added  loss,  which  is  not  a  motor  loss.   This  is  the  loss  due  to  the 
resisting  torque  of  the  water  seal  mechanism.   When  rotor  blade  pressures 
are  not  taken  this  loss  goes  on  the  data  card  as  zero. 

For  computing  purposes  the  torque  loss  can  be  converted  to  equivalent 

watts . 

• 
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W      =  Weight   on   Scales   x  Lever  Arm  x  RPS   x   2   rr     x  A  Constant 

Lever  Arm   =   6"    =   0.5' 

RPS    =  RPM/60 

2rr  -■=  6.2832 

Constant  =  1.356  (watts/f t-lbs/sec) 

W   =  Wt.  x  RPM  x  0.07095 

Ju 

4.   DATA  REDUCTION  USING  FORTRAN: 

Program  CMPRSOR  was  written  to  reduce  data  from  the  compressor 
test  rig.   A  set  of  program  cards  is  available. 

With  the  facilities  provided  by  the  school  it  is  necessary  only  to 
know  how  to  fill  out  Fortran  data  forms  so  the  data  can  be  card -punched,, 

There  are  21  items  to  be  entered  as  data.   The  data  cards  must  be 
in  the  following  order: 

First:   The  one  card  with  7  items  of  information. 

Second:   1  to  4  cards  with  "RC(J)"  values 

Third:   The  two  cards  with  the  12  data  iteiaa-  required  for  flow  and 
efficiency  computations  at  the  first  flow  rate. 

Fourth:   1  to  4  cards  with  manometer  board  readings  at  the  first  flow 
rate . 

Repeat  "Third"  and  "Fourth"  for  "N"  number  of  flow  rates. 

One  line  on  the  Fortran  form  equals  one  card,  except  that  the  Fortran 
form  has  only  12  spaces  and  the  card  has  80.   Where  all  80  spaces  on  a 
card  are  used,  two  lines  on  the  Fprtran  form  are  required.   For  this  pro> 
gram  70  spaces  (f<pr  7  items)  can  be  on  one  line  and  the  last  data  item 
can  be  entered  on  the  first  8  spaces  on  the  second  line.   Eighty  card 
spaces  are  used  nfor  entering  manometer  readings  and  the  "RC(J)"  values. 

All  data  entries  are  entered  in  10  spaces  per  data  entry  item  --  thu..' 

123456789  10 

BA  =   30.06   is  3      0.06 

RPM  =   1795   is  17      9      5. 

PZ   =  2.20846   is  2.20846 

Floating   point  numbers  may  be   placed  anywhere   in   the  assigned   space. 
They   are    shown   here    in    the    first   spaces.      T$is    system  of   palcement 
is    the  most   convenient   for    the   card  punch  operators.      A  decimal   point 
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x 

must  always   be   included   for   floating   point  numbers. 
The   data    is    coded    for   ease   of  memory   thus:: 
BA   =  barometer   reading 
EP  =   e    (epsilon) 

Six   items   bo   be   entered  on    the    first   card: 

M  =  number   of   pressure    readings    for   each    flow  rate. 

N   =  number   of    flow  rates. 

NOTE:      M  and   N  are    integers   and  must    therefore   be   entered  at    the 
right-hand   end   of    the    space    rpovided.       IF  M   =    1.0    for    10 
pressure    readings,    enter   it   as: 

123456789      10 

1        0 

There    is    no    decima|.    {ioint. 

BA  =  barometric  pressure  in  inches  of  mecury. 

PZ   =  cp0   -"-   See   p.    74  of  TR12$S  ,  ... 

u  |/„-  ^Values  depend  on  orifice 

EP  =  e  ...  See  p.  74  of  TR12  W§    diameter 

TARE  1  =  Tare  of  micromanometer  for  Pi  and  Ap  . 

TARE  2  =  Tare  of  micromanometer  for  A  42. 

Next  cards  f-n  1  to  4: 

RC(J)   Correction  for  centrifugal  force  of  air  in  tubes  for  rotor  pressures 
Coded  "Radius  Correction".   Find  values  on  p.  75  of  TR13°   Enter  "M"  values 
For  pressures  not  on  rotor  blades  enter  RC(J)  as  zero.   One  blank  card 
enters  8  zero  values.   (Blank  cards  used  here  should  not  be  confused  with 
the  one  blank  card  between  program  cards  and  data  cards.) 

For  each  flow  value  the  following  items  are  entered  on  two  successive 
cards,  the  first  7  on  the  first  card,  the  second  6  on  the  second  card. 

1.  Wl   =  Watts    reading   on    left-hand  meter 

Note:   The  meter  reading  (A,  V,  W)  are  entered  as  read. 
The  program  multiplies  by  the  correct  factor. 

2.  W2  =  Watts  on  right-hand  meter.   This  may  be  negative. 

3.  A  =  Amperes.   Enter  average  of  3  readings. 

4.  V  =  Voltmeter  reading. 
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5.  TW  =  Temperature  of  winding  f-   C. 

6.  WT  =  Weight  in  pounds  recorded  on  scales.   This  is  for  measurements  of 

rotor  pressures  only.   If  other  pressures  enter  WT  as  zero. 

7.  RPM  =  RPM 

1.  PI  S  Pj' 

2.  Tl  =  Temperature  in  discharge  pipe  --   F. 

3.  PA  =  ha   =  Atmospheric  reading  on  manometer  board  in  centimeters  of 

water . 

4.  DP  =  £p'   =  Inches  of  water. 

5.  A42  =  Plenum  pressure  reading  at  point  A42 .   Read  from  micromanometer 

in  inches  of  water. 

6.  TA  =  Temperature  of  air  entering  compressor  --   F- 

"M"  manometer  readings  are  entered  for  "N"  flow  rates .  PP(J)  on  computer 
program.   Centimeters  of  water. 


For  each  flow  rate: 

1-8   readings  -  1  card 
9-16   readings  -  2  cards 
17  -  24   readings  -  3  cards 
25  -  32   readings  -  4  cards 

The  table  of  results  is  self-explanatory.   The  pressure  (t|f)  values 
are  printed  according  to  the  order  in  which  corresponding  manometer  readings, 
PP(J) ,    were  entered  as  data. 


-9- 


TN  No.  64T-1 


PROGRAM  CMPRSOR 


DIMENSION  PSI  (96),  RC(96),  PP(96) 

Ms  N,  BA,  PZ,  EP,  TARE1S  TARE  2 


=  1 


B3SA  =*  BA*  (13.55) 

DO  10   I  -  1,  N 

READ  500,  Wl,  W2,  A,  V,  TW,  WT ,  RPM 

READ  500,  PI,  Tl,  PA,  DP,  A42 ,  TA 

TTA  -  TA+459.4 
PP1  -  Pl'TARE  1  + 

TT1  s  Tl  +  459.4 


A2  -  FZ*(1.-EP*DP/FP1) 

A3  s  1300. /RPM  *  PP1/BBA  ft  TTA/TTl 
A4  -  SQRTF(TT1/520.*DP/PP1) 
PHI  =  A2  *  A3  *  A4 

OUT  -  (A42  *•  TARE  2)  *  PHI  *  RPM 
V  =  V  *  20. 
A  »  A  *  10. 
WA  s  40.  *  (Wl  +  W2) 

A2  -  -210. +  .8929  *  V- . 0001047'.  *  V  *  V 
IF  =  (A- 12.)  1000,  1000,  2000 
1000   A3  •  2.73  *  A  *  A 


2000   A3;-  4484.-737.  *  A+33.  *  A  *  A 
3000   A3  -  A3  *  (.904  +  .00385  *  TW) 

A4  -  .000556  *  (1800. -RPM)  *  (WA^A3) 
A5.=  WT  *  RPM  *  .07095 

ETA.S  (OUT/(WA-(A2  +  A3  +  A4  +  A5)))  *  151.196 
READ  500,  (PP(J),J=1,M) 

DO  20  J  -  1,  M  „,,x 

20   PSI  (J)  s  ((PA-PP(J))  *  TTA  (BA*RPM*RPM)  *  64  10.2)  +  RC(J) 
WRITE  OUTPUT  TAPE  2,700,  PHI 
WRITE  OUTPUT  TAPE  2,800,  ETA 
WRITE  OUTPUT  TAPE  2,300 
WRITE  OUTPUT  TAPE  2,900 

WRITE  OUTPUT  TAPE  2,600,  (PSI(J),  J=  1,8) 
IF  (M-8)  10,  10,  30 
30   WRITE  OUTPUT  TAPE  2,1100 

WRITE  OUTPUT  TAPE  2,  600,  (PSI(J) ,  J-9  ,  16) 
IF  (M-16)  10,  10,  40 
40   WRITE  OUTPUT  TAPE  2,1200 

WRITE  OUTPUT  TAPE2 ,600 , (PS1(J) ,  J=  17,24) 
IF  (M-24)  10,  10,  50 
50   WRITE  OUTPUT  TAPE  2,  1300 

WRITE  OUTPUT  TAPE  2,  600,  (PSI(J),  J=  25,32) 
10   CONTINUE 

END  FILE  2 
100   FORMAT      (2110,  5F10.5) 

300   FORMAT  (27.11'.  PSI  (PRESSURE  COEFFICIENT)//) 
500   FORMAT   (8F  10.5) 
600   FORMAT   (8F9.5) 

700   FORMAT   (////26H  PHI  (FLOW  COEFFICIENT)  -  F6.5//) 
800   FORMAT   (2011  ETA  (EFFICIENCY)  =  F7.3.9H  PER  CENT//) 
9000  FORMAT   (68H    1        2        3        4       5 

1     7         J->  -10- 
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PROGRAM  CMPRSOR  (Cont'd) 


11000  FORMAT 

(69H 

1    15 

16 

12000  FORMAT 

(69K 

1    23 

24 

13000  FORMAT 

(69H 

1    31 

32 

END 

END 

25 
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SAMPLE   CALCULATION 
Horsepower   Output   of  Motor 


Measured  Test  Data 


(1)  Input  Voltage  (ave.    of   3    legs)  =  2300 

(2)  Input  i%aperes  (ave      of    3    legs)  =  10 

(3)  Input  Watts  =  31760 

(4)  Winding  Temp,  (oc)  =  20 

(5)  RPM  =  1795 

Motor  Losses  from  Curves : 

(6)  Combined  Friction  -  Windage  etc.  =   1290 
Stator  Copper  Loss  at  25°C     =  273 
Temperature  Correction  Factor   t  .981 

(9)   Stator  Copper  Loss  at  Winding  Temp.  273  x  .981  =   268 

(10)  Slip  Correction  Factor    .00834 
Input  Watts  31760 
Minus  Stator  Copper  Loss   268 

(11)  Input  Factor  31.492 

(12)  Rotor  Copper  Loss 

input  Factor  3.492  (11)  x  Slip  Corr.   .00834  (10)  =   262 

(13)  Total  Motor  Losses  (6)  +  (9)  +  (12)  1820 

Motor  Output; 

(14)  Input  31760   (3)  -  Losses  1820   (13)  =  29940 

(15)  Watts  Output  (14)  x  .001341  =  HP  =  40.2 
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Notes  on  the  USe  of  the  Compressor  Test  Rig  for  a  laboratory  test 

1.  The  oil  in  the  bearings  should  be  checked  when  starting  for  the  first 
time  after,  say,  a  week  of  no  running.   For  shorter  periods  of  tisae,  loss 

of  oil  would  be  indicated  by  pools  of  oil  underneath  the  motor.   The  bearings 
run  cool  so  any  bearing  heat  would  indicate  trouble.   The  motor  should 
take  from  three  to  five  minutes  to  come  to  a  complete  stop. 

.   There  has  been  no  bearing  trouble  of  any  kind,,   However,  there 
is  an  indication  that  bearing  friction  is  a  considerable  motor  loss.   When 
the  electrical  instrumentation  was  first  installed  it  was  found  that  the 
input  to  the  motor  decreased  during  the  first  six  hours  of  continuous  running. 
One  curve  shows  4  decrease  of  8.37o  4n  the  first  four  hours.   This  was 
analysed  as  due  to  lowered  bearing  oil  viscosity  because  of  increased  oil 
temperature.   If  possible  the  compressor  should  be  run  for  two  hours  just 
before  a  laboratory  class  ig  held. 

2.  The  motor  must  be  started  with  the  switches  to  the  low- resistance 
shunts,  that  are  across  the  instrument  leads B  closed.   If  they  are  left 

open  the  high  starting  amperage  can  ruin  an  instrument.  There  is  also  a  switch 
on  each  of  the  meters  that  must  be  closed  in  order  to  read  the  meter. 
These  switches  may  be  locked  open  or  they  may  be  held  closed  for  the  time  that 
the  meter  is  being  read. 

3.  The  sharp-edge  orifice  in  the  discharge  pipe  used  for  flow  measurements 
has  an  inside  diameter  of  13.535".   If  the  orifice  iss  as  before,  taken  out 
and  the  edge  sharpened  by  a  machining  operation ;i  the  inside  ftjameter  must  be 
measured  accurately  and  the  value  recorded  for  flow  measurement  calculations. 

4.  There  are  several  possibilities  for  pressure  leaks  in  the  lines  from  th* 
orifice  flange  to  the  micromanometer:   first,  the  connections  at  the  orifice 
flanges  may  leak  and,  second,  there  is  a  possibility  of  pressure  leaks  at 

the  quick-disconnect  arrangement  (the  brass  block  and  rubber  cork  rig)  and, 
finally,  the  connections  at,  or  even  inside,  the  micromanometer,  arap  leak. 
For  a  rough  check  on  the  possibility  of  leaks  it  is  noted  that,  at  full  flow, 
Ap  is  about  14"  of  water  for  30  conterf&ow  setting  of  the  inlet  vanes  and 
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about  10"  for  zero  setting  (vanes  radial) .   This  is  for  the  usual 
40°  setting  of  the  vanes  at  the  rotor  discharge  and  for  the  13.535"  flow 
orifice,  which  is  the  one  used.   The  13.535"  orifice  was  Atiginally 
13.500"  I.D.   The  11,5"  and  15.5"  orifices  are  not  used. 

5.  For  a  test  the  following  i Cams  must  be  at  hand 

A  thermometer 
Two  raicromanometers 
Four  clip  boards 

Blank  data  sheets,    There  are  kept  in  Room  123,  Bldg.  234. 
RPM  counter 

USNPGS  Technical  Report  No.  12 
Several  copies  of  this  publication 

Ozalid  print  of  cross-section  of  test  rig,  drawing  No^  665o 
Ozalid  print  of  cross-section  of  water  seal  mechanism,  drawing  No.  542] 
One  micromanometer  is  used  for  flow  measurement. 

One  micromanometer  is  connected  to  pressure  tap  A42  for  the  plenum 
pressure. 

A    Pressure  lines  will  be  connected  to  a  40- tube  manometer  board  as 
outlined  by  the  professor  conducting  the  class. 

The  ozalid  prints  are  to  aid  the  professor  holding  the  class  in  explaining 
the  rig  and  its  use. 

6.  When  taking  rotor  vane  pressures  great  care  must  be  taken  in  entering 
the  water  used  for  isloating  the  pressure  chambers.   If  the  water  is  let  in 
at  too  high  a  head  some  of  it  may  get  into  one  or  more  ot  the  small- 
diameter  pressure  lines  leading  to  the  tap  locations.   When  this  happens 

a  reading  at  the  tap  affected  cannot  be  taken  until  the  rotor  is  removed 
and  the  line  cleared.   The  can  that  is  raised  fo  force  water  into  the  water 
seal  mechanism  is  raiaed  Slowly  two  or  three  feet  above  the  center- line 
of  the  compressor  and  lowered  as  soon  as  all  chambers  are  isolated.   A 
manometer  board  pressure  reading  not  the  same  as  that  from  either  of  the 
two  adjacent  isolation  chamber  readings,  and  constant ,  is  an  indication 
that  sufficient  water  has  entered  to  isolate  the  chamber  connected  to  that 
particular  tube.   More  than  sufficient  water  may  cause  water  to  be  trapped  in 
the  pressure  lines  from  the  rotor  tap^lacafcions . 
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